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We calculate the longitudinal and spin-Hall conductances in four-lead bridges with Rashba - 
Dresselhaus spin-orbit interactions. Numerical results are obtained both within Landauer-Biittiker 
formalism and by the direct evaluation of the Kubo formula. The microscopic Hamiltonian is 
obtained in the tight-binding approximation in terms of the neareast-neighbor hopping integral t, 
the Rashba spin-orbit coupling Vr, the Dresselhaus spin-orbit coupling Vd and an Anderson-like, 
on-site disorder energy strength W . We reconfirm that below a critical disorder threshold, the spin- 
Hall effect is present. Further, we study the effect on the two conductivities of the Fermi energy, 
Rashba/Dresselhaus coefficient ratio, and system size. 
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I. INTRODUCTION 

Known to exist for a long timei, the spin-orbit (SO) 
coupling in two dimensional electronic systems (2DEG) 
has received a lot of attention lately motivated by its po- 
tential applications in spintronics. Recent experiments^ 
have demonstrated that the magnitude of the spin-orbit 
coupling can be modified by a voltage gate, hence gen- 
erating the premise of the possible manipulation of spin 
currents by electric fields alone. The two sources of the 
spin-orbit coupling are the inversion asymmetry of the 
confining potential in the direction perpendicular to the 
2DEG (Rashba) and the bulk asymmetry and interface 
inversion asymmetry (Dresselhaus)^. 

In a very interesting development^, Sinova et. al pre- 
dicted that a spin-Hall current of transverse spin compo- 
nent appears in a 2DEG with SO coupling as a response 
to a in-plane electric field. This spin current has a uni- 
versal value, equal to b/Stt. The intrinsic spin-Hall ef- 
fect is quite different from the extrinsic spin-Hall effecli^ 
proposed by Hirsch, which is generated by impurity scat- 
tering. The possible existence and persistence in disor- 
dered systems of the intrinsic spin-Hall effect (SHE) have 
been the focus of many recent papers^. The question of 
whether arbitrary small amounts of disorder suppress or 
not the intrinsic SHE is still awaiting a definite answer. 
Some analytical calculations^ claim that SHE does not 
survive even in the weak disorder regime, while others-: 
provide arguments that SHE is robust and weak disorder 
in the system is not enough to destroy this effect. While 
the problem was studied in more detailed using analytical 
methods, there are few unbiased numerical calculations-- 
at present. 

In this work, we present numerical results for the longi- 
tudinal and spin-Hall conductivities of a 2DEG with spin- 
orbit interactions, both Rashba and Dresselhaus, in the 
presence of disorder. These values are obtained within a 



spin-dependent Landauer-Buttiker formalism, developed 
for a microscopic Hamiltonian written in a tight-binding 
approximation that incorporates both the spin-orbit in- 
teraction and disorder. As a further check, we calculate 
the same conductances by using the Kubo formalism and 
find good agreement between the two sets of results. Our 
findings suggest that the spin-Hall effect occurs in disor- 
dered systems, for as long as the disorder remains below a 
critical threshold value. We also study the dependence of 
the conductivities on the Fermi energy, system size, and 
on the relative strengths of the two types of SO coupling. 

In the section II of the paper we present the general 
framework of spin-dependent LB formalism used for com- 
puting the spin-Hall conductance, while in section III we 
show and discuss our results. For comparison, in the ap- 
pendix, we compute the same conductances by using the 
Kubo formalism. 



II. SYSTEM DESCRIPTION 

The single particle Hamiltonian for an electron of mo- 
mentum p = {px,Py), spin <T = {ax, ay,az), and effective 
mass m*, in a 2DEG with Rashba (a) and Dresselhaus 
(/?) spin-orbit interactions is: 

H = ^— ^ -I- a {axPy - (Typx) + /3 {axPx - cryPy) ■ (1) 

The relative strengths of the Rashba and Dresselhaus 
terms, a/P describing the spin-orbit coupling in semicon- 
ductor quantum wells, are available from photocurrent 
measurementsiS. The interplay of the two SO couplings 
has been also lately subject to intense theoretical inves- 
tigations with respect to other physical phenomena such 
as magneto-oscillation phenomena in quantum wells or 
spin splitting of the electron energy states in quantum 
dotsi^. 
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We discretize the Hamiltonian using a tight-binding 
approach, where the solution domain is fiUed with a reg- 
ular virtual lattice. The Hamiltonian is constructed over 
this lattice assuming only neareast-neighbor coupling. 
This can be done straightforwardly by using the pro- 
jections on X and y direction of the momentum opera- 
tor p = —ifiSI in Eq. QJ. The resulting tight-binding 
Hamiltonian is: 



H ^ e^c\^C^o, - t ^ c\^Cja (2) 

i 

i 

Here t — /{2m*aQ) is the hopping integral, Vr — 
ha/a^) and Vd — fiP/ao are the Rashba and Dresselhaus 
coupling strengths, respectively, renormalized by the lat- 
tice constant oq, and Sx and Sy are the unit vectors along 
the X and y directions. The hopping matrix element t 
represents the unit of energy in our calculations. The 
second, third, and last terms in Eq. can be combined 
and a compact expression for the Hamiltonian can be 
written in the form: 



H 



y^£tc,;„c. 



t 



E 



(3) 



where Cj^ ("^j/^) annihilation (creation) operator of 

an electron of spin index fj, at site j. The first term in 
Eq. Q is the on-site disorder, as in the Anderson model, 
with Ei, a random energy generated by a box distribu- 
tion Ei S [—14^/2, W/2]. The SO interactions are directly 
incorporated in the hopping term which acquires posi- 
tion and spin dependence. The Hamiltonian given by 
Eq. lO is studied in a NxN square lattice, as presented 
in Figure Each metallic lead attached to the sample 
is considered a perfect semi-infinite wire, without disor- 
der and SO interactions. Vr and Vd are also assumed to 
be zero in leads 3 and 4 in order to avoid spin flips at 
the boundaries. Throughout our calculations we use the 
same values for the cross sections of leads and sample, 
in order to eliminate scattering induced by the wide-to- 
narrow geometry 

Within the LB formalism the total current in termi- 
nal p is given by Ip = e'^ /hJ^q^p'^pgi^P - ^g) where 
the sum is over all the other leads q connected to the 
system. Spin current can be defined in a similar way, 
up to a constant: J«p» = e/ (47r) 7;^<r(^P - Vq). 

The voltages are computed by considering ballistic trans- 
port between all the connected terminals and impos- 
ing the following boundary conditions: V2 = (fixes 
the arbitrary zero of voltage), = "^e/hj^a -^s^a" — 0, 
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FIG. 1: Graphical depiction of the lattice model used for 
computing the spin-Hall conductance. Four metallic leads 
(represented as the dashed regions) acting as injector (1), de- 
tector (2) and voltage probes (3 and 4) are are attached to 
the 2DEG. Position and spin dependence are not explicitly 
decided for the hopping integral. 



I4 = 2e/?i^^ J^^" = (terminals 3 and 4 are voltage 
probes) and Ii + I2 = 0, (guarantees that current fiows 
between terminals 1 and 2). The zero temperature con- 
ductance, G, that describes the spin-resolved transport 
measurements, is related with the transmission matrix 
T, as in : 



G 



h 



(4) 



[Indices p and q were suppressed in writing Eq. |0J] . T^^" 
represents the transmission probability over all the con- 
duction channels to detect a spin /x in the lead p aris- 
ing from an injected spin v electron in lead when 
both spin-fiip and non-spin-flip processes are consid- 
ered. The transmission coefficient can be calculated as 

= rr[r^Gfir^GA] where T^ = ^(E^ - E^t) ^ith 
the retarded self-energy due to the interaction between 
the sample and the lead for spin-channel /i. The self- 
energy contribution is computed by modeling each termi- 
nal as a semi-infinite perfect wire. In our tight-binding 
model, the hopping between the lead orbitals and be- 
tween the leads and the sample orbitals are equaU^ with 
t (unit of energy). The self-energy matrix, which is diag- 
onal in spin indices, can be written as: 



E - ( 

^" " ' eJ, 



(5) 



with Ej, — Ej, for a perfect metallic lead. The re- 
tarded Green's function is computed as Gr = {Ep — H — 
X]p=i where Ep is the Fermi energy and H is the 

Hamiltonian in Eq. (|3Jl. The advanced Green's function 
is, of course, Ga = G\. 

In the LB formalism, the total scattering between two 
lead p and q can be simply written as the sum over all 
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spin components T, 



pq 



useful combinationsi^ are TJ 



and Tp^* 



— TTT -I- T-LT 

pq ' pq 



Two other 



pq ' pq pq PQ 

pg pq ' pq PQ 

Tjj - Tpi^l. r°«* represents 
the difference between the transmission probabiHties to 
detect an electron in the lead p arising from an injected 
spin f (j) electron in lead q. These expressions allow us 
to compute the spin-Hall conductance, as 



H 



j-sptn 
^3,T 



jspin 
^3,1 



Vi 



(6) 



Finally, by using the voltages obtained inverting the mul- 
tiprobe equations, the spin-Hall conductance becomes: 



G. 



H 



e/(87r)(r, 



out 
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TiOUt 
43 



rpOUt 

^23 



r--2T37). (7) 



At the same time, the longitudinal conductance, Gl — 
h/iVi — V2), is written as: 



Gl = e^h (T21 + 0.5 T32 + 0.5 T42) 



(8) 



when four terminals are connected to the sample as in 
Figure ^ The spin-Hall and longitudinal conductances 
are the central quantities of our analysis. In the next 
section, we present results showing their dependence on 
the Fermi energy, system size, and disorder strength. 



III. RESULTS AND DISCUSSION 

Clean limit. The clean limit dependence of the spin- 
Hall conductance (SHC) on the Fermi energy is shown 
in Figure HJ The electron-hole symmetry is preserved 
throughout the calculation, so the SHC vanishes at the 
band center Ep — and is an odd-function relative to 
the Fermi energy, in agreement with the results of Ref . 
The small oscillations observed in the energy dependence 
are finite size effects related with the discontinuities in 
the self-energy contribution from the terminals and with 
the discrete energy levels. 

Another important parameter is the ratio r = Vr/Vd- 
When 7- = 1, for any energy, GsH = (see Fig. 12 
- right panel, and Fig. O. For a hole- like behavior 
{Ep < 0) and r > 1, GsH is positive, while for r < 1, 
GsH changes sign, demonstrating that the spin current 
is generated in the direction of the major driving fielcii^. 
Experimentalljiifliii, the tuning parameter r could be 
varied between 1.5 and 2.5. 

Figure |2| presents the effect of the Dresselhaus SO 
coupling on the longitudinal conductance as function of 
Fermi energy, for a fixed value of the Rashba coupling. In 
contrast, in Fig. 2]the Fermi energy is fixed to Ep — ~2t 
and the longitudinal conductance is plotted as function of 
both Rashba and Dresselhaus interactions. Here, r — 1 
still represents a symmetry line in the parameter space 
(Vr, Vd)- For a fixed value of the Fermi energy, we found 
the symmetry relation: Gl{Vr, Vd) = Gl{Vd,Vr). 

In Figure |S1 we present the spin-Hall conductance as 
function of Vr and Vd ■ SHC is anti-symmetric along the 
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FIG. 2: Upper left panel: The Fermi energy dependence of 
the spin-Hall conductance (SHC) of a two dimensional four- 
probe bridge in the clean limit, for a fixed Rashba coupling 
Vr = 0.5 and for different Dresselhauss energies. Upper right 
panel: SHC dependence on Yd for different Rashba couplings, 
for Fermi energy Ep — —2t in the clean limit. For Vr = Vd, 
the SHC vanishes. The system size is 20 x 20. Lower panel: 
The SHC represented as function of Ep for Vr = 0.06 and 
Vd = 0.0. 




FIG. 3: The longitudinal conductance as function of the Fermi 
energy for different Dresselhaus SO couplings, as indicated. 
The system size is 20 x 20. 



Vd = Vr line. The Fermi energy is fixed at Ep — ~2t 
and the system size is 20 x 20. For a lattice parameter 
of oo = 5.0 nm and electron effective mass m* = 0.068m 
(in GaAs), the hopping integral is i ~ 19.0meV. A typi- 
cal value for the Rashba coupling^^ is ~ 50 — 80 meV-A, 
which corresponds to Vr = 1 — 1.6 meV, with a typical 
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FIG. 4: (Color online) Longitudinal conductance plotted as 
a function of Vr and Vd for a system size 20 x 20 and for 
a Fermi energy Ef = —2t. The spectrum is anti-symmetric 
along the Vd = Vr line. The spin-Hall conductance is positive 
for Vr > Vd, negative for Vr < Vd and vanishes for Vd ~ Vr. 
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FIG. 5; (Color online) The spin-Hall conductance plotted as 
a function of Vr and Vd for a system size 20 x 20 and for 
a Fermi energy Ef = —2t. The spectrum is anti-symmetric 
along the Vd = Vr line. The spin-Hall conductance is positive 
for Vr > Vd, negative for Vr < Vd and vanishes for Vd = Vr. 



ration Vfj/t ~ 0.05 — 0.08. The results presented in Fig. 
121 (upper panel) are beyond the experimental reach. In 
Figure 121 (lower panel) we represent the Fermi energy de- 
pendence of the SHC with a experimental accessible value 
for the spin-orbit interaction strength, and, as expected, 
the SHC amplitude is strongly reduced. However, the 
overall behavior is preserved. 

The effect of scaling as function of system size is pre- 
sented in Figure Spin Hall conductance is essentially 
constant up to, at least, system sizes 50 x 50. However 
we emphasize that the effect of boundaries, due to the 
attached leads, may be very important and in principle 



can hide the true nature of the bulk spin-Hall effect. 

Our analysis shows that in the clean limit, a non- 
universal value for SHC exists, in agreement with other 
numerical calculations^. SHC strongly depends on the 
strength of spin-orbit couplings, while the spin current 
is always along the driving field in the system and de- 
pends on the relative strength of the Rashba and Dres- 
selhaus couplings. In the hole-like (electron-like) regime, 
characterized by Ep < {Ep > 0), sgn{SHC) ~ 
±sgn{Vi, - Yd). 

Arbitrary disorder. A system with time reversal 
symmetry, but with spin rotational symmetry broken by 
the spin-orbit coupling, belongs to the symplectic uni- 
versality class. It is well established by now that SU (2) 
models with chiral symmetry exhibit an Anderson tran- 
sition in two-dimension^^. Critical disorder strength was 
estimated to be Wc — 5.9 and the critical exponent for 
the localization length i' ~ 2.74. In our model, different 
values for the hopping coupling may lead to different val- 
ues for the disorder strength. However, it is understood 
that SHC cannot survive in the insulating regime of a 
2DEG, because any localized state cannot contribute to 
SHC. It is still not clear whether SHC vanishes in the 
diffusive transport when the mobility edge ±Ec moves 
towards the band center and localized states in the band 
tails coexist with extended states in the band center. To 
answer this question we study the effect of disorder on 
SHC. In Figure [3 (Left panel) we represent the SHC as 
function of disorder strength for different Yd- We find 
that GsH can be suppressed by a strong scattering when 
W > 4 — 5, close to the metal-insulator transition disor- 
der strength. In the left panel the Dresselhaus coupling 
is zero and the Rashba coupling strength dependence of 
SHC is presented. For comparison we have plotted also 
the result when no disorder is present in the system. 

Energy dependence was also considered in the presence 
of disorder, (see Figure [7| lower panel). 

In the insulating regime, all states are localized, so 
the absence of extended states available for transport at 
the Fermi level leads to a vanishing SHC. When disor- 
der is weak enough, extended states in the band cen- 
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FIG. 7: (Color online) Left: Disorder strength dependence 
of the SHC with Rashba spin-orbit strength Vr = 0.4. Right: 
SHC as function of Rashba coupling for different disorder 
strengths W. Lower panel: SHC as function of Fermi en- 
ergy and for different disorder strengths W. Electron-hole 
symmetry is preserved in the presence of disorder. In the up- 
per panel Fermi energy is Ep = —2.0t. System size is 16 x 16 
in all figures. 

ter coexist with insulating states localized mostly in the 
band tails. These extended states may be responsible for 
non-vanishing SHC when small amounts of disorder are 
present in the system. 

It is well known that in Landauer-Biittiker formalism 
the attached leads play an important role, affecting the 
system self-energy, and can alter the nature of the bulk 
spin-Hall effect, while this is not the case in the Kubo 
formalism. 

To study the effect of terminals on the spin-Hall and 
longitudinal conductances we did a direct calculation of 
conductances within the Kubo formalism (see the ap- 
pendix for further details). We found good agreement 
between the conductance values obtained in both the LB 
and Kubo formalisms. In Figure ^we present the these 
results for a system of size 16 x 16. The electron-hole 
symmetry is also preserved in the Kubo formalism, so 
the spin-Hall conductance vanishes at half filling, as in 
the LB formalism. 

In Ref. 16, Sinitsyn et. al. and Shen use the Kubo 
formula to compute the spin-Hall conductance analyti- 
cally, when both the Rashba and Dresselhaus couplings 
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FIG. 8: (Color online) Longitudinal conductance as function 
of disorder strength (upper left panel), and as function of 
Rashba spin-orbit coupling strength (upper right panel) for 
different disorder amplitude. Lower panel: Effect of disorder 
on the energy dependence of the longitudinal conductance. 
The system size is 16 x 16. Disorder average is over 1000 
samples. 

are considered. As in our case, they find that the spin- 
Hall conductance vanishes when the Rashba and Dres- 
selhaus couplings have the same strength. The predicted 
value of the SHC, however, is a constant zLe/Stt, depend- 
ing on the ration Vr/Vd- In contrast, in our numeri- 
cal approach, the SHC is no longer a universal constant, 
but rather a function of the Fermi energy and of the 
Rashba/Dresselhaus coupling strengths. 



IV. CONCLUSIONS 

In this work we have investigated the longitudinal and 
spin-Hall conductances of a two dimensional electronic 
system with Rashba and Dresselhaus spin-orbit coupling 
in the framework of a tight binding approximation. For 
the main part of the work we have used Landauer- 
Biittiker formalism combined with Green's function ap- 
proach to study the efect of spin-orbit coupling and disor- 
der on Gi and GsH- Our results for the spin Hall conduc- 
tance, as function of Fermi energy and disorder strength, 
in the case when Dresselhaus coupling is neglected, agree 
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with the results of Ref. '5 which is a special case of the 
present model. 

We have also computed the Fermi energy dependence 
of longitudinal and spin-Hall conductances in the Kubo 
formalism. The good agreement found between the two 
sets of conductances computed in LB and Kubo for- 
malisms strengthens the assumption that the spin-Hall 
effect is a bulk property of the system. However, further 
studies are needed in order to clarify the role of terminals. 
For example, one can investigate the scaling of spin-Hall 
conductance as function of the system size, both in the 
Landauer-Biittiker and the Kubo formalisms. 
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APPENDIX A: COMPARISON WITH THE KUBO 
FORMALISM 

In this Appendix we present the derivation for the 
Kubo formula used for computing the longitudinal and 
spin-Hall conductances. 
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FIG. 9: (Color online) Longitudinal (upper panel) and spin- 
Hall conductance (lower panel) for a system size of 16 x 16. 
Solid lines represent results obtained using Kubo formalism 
while the dashed lines are obtained using Landauer-Biittiker 
method. Results are average over 500 samples. Energy is 
measured in units of t. 

In terms of single particle states, the longitudinal con- 
ductivity can be computed from the general Kubo for- 



(7L (r,r') ^ ~ihY^ 



Pn' - Pn jn' \jxir)\ n) {n \va;{r')\ n') 
En, - En E,,, -En+i m+ 

(Al) 

Similarly, for the spin-Hall conductance we write: 

Pn' - Pn Im {n' |j|(r)| n) {n \vyir')\ n') 



asH (r,r') = ^ 



En 



En 



En' — En 



(A2) 

The single-particle states are constructed from the site 
orbitals as 5]^ = J2i a '/'" ^) ^L- Operator b}^ stands for 
the creation of a single particle state \n) from the one- 
electron wave functions, V'n(*jQ^)- The wave functions 
Tpn ih ck) and the corresponding eigen-energies En can be 
easily obtained by solving the eigen- value problem for the 
Hamiltonian J^J. 

The velocity operator is defined by the commutator: 
ih'v — [r, H], while the spin current is given in terms of 
the anticommutator between the velocity operator and 
Pauli matrix az- j| = fi./4 {cr^, Wj.}. A simple quantum 
mechanics calculation gives for the current and for the 
spin-current operators the following expressions: 



(n |v| n') 



1 

i h 



a, 13 



(A3) 



(A4) 

In Eq. (EH, H = {az(^ l,H}. 

At T = OK, when the Fermi function derivative is 
approximated by a delta function, we write: 



Pn 



Pn 



En 



En 



-&{En - Ef). 



E) 



(A5) 



Incorporating Eqs. (|X3|) . HA4|1 . ((X5|) in Eqs. ((XTI) and 
(IA2|) a simple expression for longitudinal and spin-Hall 
conductance in terms of single-particle wave functions 
and eigen-energies is obtained. We note that only terms 
at the Fermi levels give contributions to the longitudinal 
conductance therefore we keep only the delta function 
part from {En' - En + i ?iO+)~^, in Eq. HAll) . In this re- 
spect, the longitudinal conductance is a sum of weighted 
delta functions which have to be broadened into func- 
tions having a finite width (for example a Lorentzian). 
When the spin-Hall conductance is computed, the prin- 
cipal value of {En' - En + i ?iO+)~^ is needed in Eq. HA2|I. 
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